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Abstract One central challenge in high-T c superconductivity (SC) is to derive a detailed 
understanding for the specific role of the Cu-d x i_ y i and 0-p X]V orbital degrees of freedom. 
In most theoretical studies an effective one-band Hubbard (1BH) or t-J model has been 
used. Here, the physics is that of doping into a Mott-insulator, whereas the actual high-T c 
O cuprates are doped charge-transfer insulators. To shed light on the related question, where the 

material-dependent physics enters, we compare the competing magnetic and superconducting 
phases in the ground state, the single- and two-particle excitations and, in particular, the 
pairing interaction and its dynamics in the three-band Hubbard (3BH) and lBH-models. 
Using a cluster embedding scheme, i.e. the variational cluster approach (VCA), we find which 
frequencies are relevant for pairing in the two models as a function of interaction strength and 
doping: in the 3BH-models the interaction in the low- to optimal-doping regime is dominated 
by retarded pairing due to low-energy spin fluctuations with surprisingly little influence of 
inter-band (p-d charge) fluctuations. On the other hand, in the lBH-model, in addition a 
part comes from "high-energy" excited states (Hubbard band), which may be identified with 
a non-retarded contribution. We find these differences between a charge-transfer and a Mott 
insulator to be renormalized away for the ground-state phase diagram of the 3BH- and 1BH- 

i i models, which are in close overall agreement, i.e. are "universal". On the other hand, we 

expect the differences - and thus, the material dependence to show up in the "non-universal" 
finite-T phase diagram (T c -values). 
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Many aspects of the physics of high-temperature superconductors (HTSC) remain mysterious, despite 
impressive progress both on the experimental and theoretical front. One key issue, is why are the HTSC 
materials composed out of CuO-i planes and what is the specific role of Cu-d and O-p orbital degrees of 
freedom? The cuprate materials in the undoped, i.e. "half-filled" situation, are "charge-transfer" insulators 
[1 2 . This fact induces an experimentally observed asymmetry between hole (h)- and electron (e)-doping: 
while doped holes go onto O orbitals and may be bound to Cu spins to form "Zhang-Rice" singlets, 
doped electrons reside mainly on the Cu orbitals |3I4I5I6I7I8I9I10| . This is believed to be intimately 
related to the more extended stability of antiferromagnetic (AF) behavior as a function of electron 
doping compared to that of hole doping. While introducing electrons on the Cu sites merely produces 
a dilution of spins, holes on O sites create a ferromagnetic coupling between neighboring Cu orbitals, 
which is significantly more effective in destroying the AF order. Beyond magnetism, the asymmetry is 
also exposed in the superconducting (SC) behavior with the /i-doped materials exhibiting SC usually 
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over wide doping regimes and with high T c 's up to 150K, whereas in the e-doped system T c is low and 
confined to a very narrow doping regime. 

In much of the high-T c theoretical studies the starting point has not been the above three-band Hub- 
bard (3BH) but instead the one-band Hubbard (1BH) and t-J models |11I3| . Here, the oxygen degrees 
of freedom are eliminated approximatively via the Zhang-Rice construction |3j. The mapping has as a 
crucial consequence that the undoped model systems are Mott insulators and no longer governed by 
the charge-transfer energy A p d between the Cu and O orbitals. Analytical and, in particular, numerical 
calculations based on these two-dimensional (2D) single-band models have demonstrated Fermi surfaces, 
single-particle spectral weights, AF spin correlations and d x i_ y i pairing correlations in qualitative agree- 
ment with experimental measurements [12113114115116117118119120] . This fact has significantly contributed 
to the wide-spread belief that the physics of HTSC is that of "doping into a Mott insulator" [2T| . However, 
how can this picture be reconciled with the charge-transfer insulator picture embedded in the 3BH-model? 
Can it be that the charge-transfer energy A p d in the 3BH-model, the size of which is already decisive for 
the accuracy of the 1-band reduction |3] , plays the role of an effective on-site Hubbard U in the 1-band 
models? 

Quantifying these ideas requires solving the strongly correlated electron problem for the 3BH- and 
1BH- models at very low energy (and/or temperatures). Early Quantum-Monte-Carlo (QMC) calculations 
for the 3BH-model showed, that characteristic features such as the doping dependence of the electronic 
single-particle excitations and their interplay with magnetic excitations are in accord with experiment 
|4I5| . However, the very low T- or ground-state properties including the SC state, could not reliably 
be resolved, due to the well-known "minus-sign" problem [22J. Embedded cluster techniques provide a 
controlled way to approach the infinite-size (and, thereby, low-energy) limit. Recently, the variational 
cluster approach (VCA) which was proposed by M. Potthoff and our group |23I24| has been shown for 
the 1BH model to correctly reproduce salient features of the ground-state (T = 0) phase diagram of the 
high-T c cuprates |14I16I17I25] , In particular, the AF and d-wave SC ground-states were found in doping 
ranges qualitatively in accord with experimental data for both e- and /i-doping such as the different 
stability of the AF phase. It can be accounted for by a "simple" 1BH model, in which the e-h symmetry 
is broken by a longer-ranged (next-nearest) hopping term. 

Regarding the possibility of the reduction to a one-band model, there are questions concerning the 
direct applicability of the Zhang-Rice (ZR) construction. As discussed in the literature before |2], the 
natural tendency of a finite oxygen band width is to delocalize and to destabilize the ZR singlets. Secondly, 
the pragmatic finding that a t-t'-U lBH-model with a significant value of t' , captures basic physics of 
the cuprates and in particular their e-h asymmetry, cannot be accounted for in a strict ZR picture [26 . 
In this picture, next-nearest neighbor hoppings are very small compared to nearest-neighbor terms (if 
again oxygen-oxygen hopping t pp is taken into account). What we will show in our calculations is that the 
3BH-model and a single-band t-t'-U Hubbard model with a significant value for t' (which may be taken 
as an empirical parameter adjusted to fit the Fermi-surface topology |27| ) exhibit a similar low-energy 
(here specifically, T = 0) physics concerning the qualitative behavior of the ground-state phase diagram 
[28J and of the single-particle excitations including the asymmetry as a function of h- and e- doping. 

Then, where do the cuprate material-specific properties appear? To shed light on this question, we 
start out in section III with contrasting the low-T c pairing mechanism with a discussion of what we know 
about the high-T c mechanism. In particular, we take up an issue raised recently by P.W. Anderson |29| . 
whether the pairing interaction in the cuprate HTSC should be considered as arising from a "pairing 
glue". Anderson argued that there is no room left for a pairing glue: while the low-T c SC pairing is due to 
the well-known dynamic screening mechanism which contains a (phononic) "pairing glue", in the high-T c 
cuprates and other unconventional SC another mechanism is at work. It may be termed "anisotropic real- 
or momentum-space" pairing. Here, the electrons pair in an anisotropic wave function (such as d-wave), 
which vanishes at the repulsive core of the Coulomb interaction. While this anisotropy is certainly realized 
for the cuprate HTSC and embedded in the momentum dependence of the d-wave gap function, there 
is, additionally, a frequency dependence of the SC gap, which tells us about the dynamics of the pairing 
interaction [30 . We study in section IV the dynamics of the SC gap function and the building up of a 
"pairing glue" which contributes to the formation of Cooper pairs in the 3BH-model and contrast it with 
the corresponding "pairing glue" for the lBH-model. The "pairing glue" question is a natural one to ask if 
one has two-particle excitations, which are responsible for the dynamics of the pairing interaction, whose 
energy scale is small and of the order of the SC gap energy. For the 1-band models this "pairing glue" 
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issue has recently been investigated [30 31 32 . Spin fluctuations create a retarded interaction, which is 
dominating the low-energy pairing. However, also a pairing contribution due to interband (lower and 
upper Hubbard band) transitions was identified [30 in the 1BH- model case. In our calculation for the 
lBH-model, this "non-retarded", i.e. spin-fluctuation dominated part is comparable in magnitude to a 
second part: For large U's this creates an "instantaneous" pairing interaction eventually going over into 
the instantaneous part previously identified in the t-J model |30| . 

There are even qualitative differences in the dynamics of the pairing interaction in the 3BH-model: 
here, the dynamics of the pairing interaction is dominated by the retarded contribution due to spin 
fluctuations. Higher-energy p-d charge fluctuations give rise to a surprisingly small further (retarded) 
contribution. As expected, this charge-fluctuation contribution gains slightly more weight, when going to 
higher hole dopings. In short, one may term the dynamics of the pairing interaction in the 3BH-model 
and, thus, in a doped charge-transfer insulator as retarded and being due to a low-energy pairing glue, 
whereas in the lBH-model and, thus, a doped Mott insulator both retarded and non-retarded parts 
contribute. The differences in the dynamics of the pairing interaction, which are due to the material 
dependent CuCVphysics, appear at relatively large energies (of order A p d) which are significantly larger 
than the SC gap. In the ground-state phase diagram these "high-energy" differences are renormalized 
away (see, Fig. [5J in section IV). We do, however, expect that they play a role in the finite-T phase 
diagram (T c - values). 

Also the low-energy single-particle excitations of the 3BH- and lBH-model are qualitatively similar, 
when comparing the 3BH-model with the (as discussed above) empirical t-t'-U model. A detailed discus- 
sion of this was presented in Ref. [25 . This concerns, in particular, where holes or electrons, respectively, 
go when one dopes away from half-filling. The corresponding nodal and anti-nodal doping regimes are 
directly related to the observed asymmetry in the robustness of the AF order. Closer inspection, however, 
also reveals differences for higher energies of 0(t p d). An example is the recently much discussed "waterfall" 
structure or high-energy kink appearing as an abrupt change in the band dispersion, which falls vertically 
at binding energies below sw OAeV [33]. This "waterfall" behavior is found to be rather pronounced in the 
3BH model, but not in the lBH-model. 

Finally, in section V, we describe a theory we have recently developed of two-particle (2-p) excitations 
and, in particular, the magnetic susceptibility in high-T c cuprate superconductors [34 . As just discussed, 
2-p excitations such as spin and charge excitations, play a key role in the pairing theory. There, they 
are indirectly embedded in the anomalous part of the self energy and the dynamics of the gap function. 
However, for an unambiguous identification, for example of a "pairing glue", one obviously needs a di- 
rect scheme for calculating these 2-p excitations in the experimentally relevant strong-correlation regime. 
Therefore, we have recently extended the Variational Cluster Approach to these excitations [31]. We found 
for the lBH-model, that the corresponding magnetic susceptibility \s reproduces salient features of the 
experiment such as the celebrated neutron-scattering resonance in the hole-doped case. Previous descrip- 
tions of the magnetic resonance have been obtained by weak-coupling [35 and/or semiphenomenological 
approaches [36 and, quite independently, by the SO(5)-symmetry argument |37| . The infinite-lattice limit 
is crucial to obtain the magnetic resonance which may be considered as a "fingerprint" of the AF order 
in the SC state. Only then are we able to differentiate between the competing AF and SC orders in the 
phase diagram. Therefore, this limit, which is obtained with the VCA, has also to be embedded in a 
controlled description of the corresponding susceptibilities. 

In section V, we also study the emergence of a similar magnetic resonance in the SC state of the 
electron-doped cuprate superconductors. We find again for the lBH-model case, that the experimen- 
tally observed resonance peak is consistent with an overdamped magnetic (spin S=l) exciton located 
near the particle-hole continuum. Because of computational reasons (reference cluster size) we can only 
demonstrate the magnetic resonance for the lBH-model. 



II) The Models 



We consider the following version of the three-band Hubbard model: 
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where d- creates a hole with spin a in the Cu-d x 2_ y 2 orbital at site i with occupation number nf a . 
Correspondingly, p'- creates an 0-p x _ y hole at site j with occupation number n^ CT . £y stands for the p-d 
hybridization tc u ~o & n d ijf for the direct O-O hopping, where the orbital phase factors are included 
[4 5 . The local orbital levels are given by e p and e^ and the charge-transfer energy is A — e p — td- Udd an d 
U pp are the Hubbard couplings on the Cu and O sites, respectively. Finally, U p d is a repulsive interaction 
for holes on Cu and O. We take typical parameters consistent with earlier extensive Quantum-Montc- 
Carlo (QMC) |4I5| as well as cluster [57] calculations. In units of the Cu-O hopping t p d = 1, A = 3, 
Udd = 8, Upd = 0.5, Upp = 3.5 and t pp = 0.5. 

The lBH-model is defined as usual, i.e. 



H= J2 t ljC \ a c ]a + uYn Anii . (2) 

a,<ij> i 

Here Cj a and c^ are annihilation and creation operators and ty denote the nearest (i„„ = 1, energy 
unit) and next-nearest neighbor (t nnn = —0.3) hopping matrix elements and U (U = 8) the on-site 
Hubbard repulsion. 

An important remark for the comparison of the two models is that while for the 3BH-model the unit 
is tpd = 1 and t nn = t = 1 for the lBH-model, when going back to eV's one has to set t p d — 1.5[eV] and 
t = 0A[eV}. 



Ill) Low-T c versus high-T c superconductivity and the question of the pairing glue. 

In a recent article [2S], P. W. Anderson (PWA) has argued that pairing in conventional "low-TV super- 
conductors (SC) has a rather different microscopic origin from that in high-T c cuprates and many other 
unconventional SC. In either case, the paired electrons have to avoid the strongly repulsive bare Coulomb 
interaction. 

In a low-T c SC this repulsion can be eliminated in favor of electron-pair binding via "dynamic screen- 
ing", i. e. 

q 2 e{(\,u)) 

where V(q, w) is the Fourier transform of the el-el interaction in both space and time and e(q, ui) is the 
dynamic screening due to both ions and electrons, i.e. e(q,u>) — £i on (q, lo) + e e ;(q, ui). 

PWA then gives an elegant discussion why, in his opinion, in the high-T c SC another pairing mechanism 
is at work, which may be termed "anisotropic momentum (or real-space) mechanism": here, the strongly 
repulsive (short-range) part of the Coulomb interaction is avoided by a mechanism suggested by Pitaevskii 
[38J and Brueckner et. al. [39 of choosing the pair state orthogonal to the repulsive core of the Coulomb 
interaction, i.e. putting the electron pairs in an anisotropic wave function (such as d-wave), which vanishes 
at the core of the Coulomb interaction. 

While we follow the arguments of PWA for the low-T c SC, we would like to present here evidence, 
that in the high-T c materials both the anisotropic "momentum space" mechanism and the dynamics are 
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Figure 1. The effective pairing interaction in the weak-coupling low-T c situation compared with high-T c QMC 
simulations for the 2D f-band Hubbard model by the Scalapino group [43 . 



at work. In what follows this will be shown for both of the above Hubbard-type of models (see chapter 
IV). This gives us, in particular, a possibility to shed light on the question to what extent both models 
lead to a similar pairing mechanism. For the 1-band Hubbard model, the reader is also referred to recent 
studies by Maier, Poilblanc and Scalapino [3U] as well as by Kyung, Senechal and Tremblay [3T|, where 
also the pairing dynamics has been studied. The above questions are clearly of relevance for the main 
issue raised by PWA, namely "is there glue in cuprate superconductors?", which is basically a question 
about the dynamics of the pairing interaction. As argued by Maier et al., if the dynamics of the pairing 
interaction arises from virtual states, whose energies correspond to the "high-energy" Mott gap, and give 
rise to the exchange coupling J, the interaction is instantaneous on the relative time scales of interest. 
However, if the energies correspond to typical "low-energy" spin-fluctuation (or phonon) excitations, then 
the interaction is retarded. In this case it makes sense to use the terminology "spin-fluctuation glue", 
which mediates the d-wave pairing. Here, we present results from numerical (variational cluster) studies, 
which provide insight into this question. 



For this presentation it is useful to step by step contrast the construction of the effective el-el interac- 
tion for a weak- coupling low-T c SC with the effective interaction of a high-T c (i.e. Hubbard type) SC. In 
the "dynamic screening" mechanism of eq. (3) one can safely replace the electronic screening by a static 
one (i.e. £ e z(q, w = 0)). This is due to the fact that typical el-energies are of "high energy" (~ 0{eF erm i) 
compared to kgT c , i.e. the SC energy scale: the plasma of other electrons then damps away the long- 
range (l/r)-behavior and leaves a Thomas- Fermi screened core e e~ xr /r (x: Thomas-Fermi constant). 
This gives rise to an essentially instantaneous interaction, which is still "very repulsive", and which - when 
averaged over the Fermi surface - is termed /j, (see below). On the other hand, for the phonon case, the 
screening acts anti-adiabatically, i.e. £i on (q, w) is dynamic, since typical phonon frequencies are of the 

2 

order of knT r . In other words, the final Fourier-transformed effective interaction is I4ff = tt^t-, — §w — vr< 
and the effective electronic interaction is screened (anti-adiabatically) by the phonon polarizations. 

Fig. [T] summarizes this weak-coupling low-T c situation and rewrites the effective pairing interaction as 
the usual sum of the dynamic phonon exchange and the Thomas- Fermi screened Coulomb part. Here uj q 
denotes phonon (ph) frequencies, <? q the el-ph coupling and uj m are Matsubara frequencies. When the two 
terms in the sum are averaged over the Fermi surface (the brackets in Fig. ^ for V g denote an average 
of the momentum transfer q over ep ), then one finds that the corresponding averaged V e g is larger than 
zero, i.e. -(A — //) > and, thus, is still repulsive. The net interaction is thus repulsive even in the phonon 
case, a fact which is required to guarantee the stability of the solid. 
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Figure 2. left-hand part: Dynamics of the pairing interaction for the low-T c case: the emergence of the difference 
in frequency scales, i.e. slow ionic-lattice polarization and fast Coulomb repulsion. Right-hand part: the SC gap 
function and its dynamics (both Re- and Im- parts) for the 3BH-model. Note that the Re-part displays a certain 
similarity to the gap function shown for the low-T c case in the left-hand part of Fig. [3] This is true for the spin- 
fluctuation contribution appearing around the characteristic energy coo + Ao = $ (see also Fig. pi)). At higher 
energies of order w — 1 and above additional structure appears mainly due to p — d charge fluctuations (see again 
Fig.|6f)). 



So, how can one ever arrive at bound pairs, if the interaction is never attractive? As is well-known 
since Cooper's seminal paper [30] preceding BCS-theory, this is due to the difference in frequency scales 
or, equivalently, in time scales of the two parts of the interaction. This is pictured in Fig. [2j the "first" 
electron (anti-adiabatically) polarizes the ionic lattice and sets up a net negative charge polarization in 
its vicinity. This first process is "slow" and happens on the frequency scale ui q of the ionic vibrations. A 
"second" electron feels this polarization, but can only profit from it and build up an effective attraction 
when the first electron "instantaneously" moves so as to avoid each other. Thus, the "high-energy" part 
of the Coulomb interaction acts only over a short time, given here by a (5-function, while the attractive 
el-ph interaction is retarded by the much slower lattice response. In other words, if the electrons forming 
the pair correlate themselves in time to avoid the short-time Coulomb repulsion, they can take advantage 
of the attractive el-ph mediated interaction and form a Cooper pair. As is well-known, this kind of 
renormalization, i. e. integrating out the "high-energy" degrees of freedom in the relative pair wave 
function (so that one arrives at one and the same low-energy cutoff of order of the Debye frequency ujd 
in Fig. [2|, can be done a variety of ways: by a ladder sum renormalization or a "pseudozation" of the 
interaction , replacing [i by \i* |41) . Now, we can have — (A — fi*) < 0, i.e. a situation where V e ff can be 
attractive, leading in Cooper's sense to the pairing instability at very low energies. 



We are interested in the dynamics of the pairing interaction in the high-T c cuprates and, in particular, 
how there the two electrons about to form a pair can avoid each other - and thus weaken their repulsion 
- by modifying the high-energy parts of their relative wave function. For this, it is instructive to look 
first at the frequency dependence of the gap function of a low-T c SC found by solving the Eliashberg 
equations [32], following a discussion by D. J. Scalapino [33] (see the low-T c case in Fig. [3]): The real part 
of the corresponding gap function A(uj) first increases as the typical phonon energy w q is approached. 
At this characteristic "glue energy", it changes sign and remains negative out to very large values of 
to. This latter observation corresponds to the instantaneous part. It just reflects the fact that the two 
electrons making up the pair avoid "short-time close-range encounters". This can be summarized (see 
Fig. p| in a kind of orthogonality relation, with the pair wave function A(oj) being orthogonal to the 
"core" (i.e. the Thomas- Fermi screened short-range part) of the Coulomb interaction [43 . The essence of 
this orthogonality relation is that in practice (47re 2 /(g 2 +x 2 )) can be replaced over the frequency integral 
A to several times the Debye frequency curj by the weak pseudopotential \x* |3T] |43| . 
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Figure 3. Avoiding the strongly repulsive part of the Coulomb repulsion: Analogy of the gap function A in the 
low-T c (w — space) and high-T c (momentum space) situation. The figure for the momentum-space gap function in 
the high-Tc- case is again taken from [43] . Here, the shaded region is the normal-state Fermi surface and the solid 
curve surrounding the Fermi surface gives the size of the d a; 2_„2-gap Ak at the momentum k on this surface. 



So then, how do the electrons in high-T c materials act so that they are seldom or never in the same 
place at the same time? 

Here, it is useful to first look at the effective interaction V e g in real-space (or, more precisely, at 
the real-space Fourier transform of the singlet vertex r(l,u> m — 0) [42 versus the separation 1 between 
the electrons in pairs). Figs. [l]&[3]plot the results for V e s obtained from QMC simulations for the 2D 
1-band Hubbard model by the Scalapino group |H| . We have obtained a similar pattern for the effective 
interaction in the 3BH-model employing also QMC, when again the two electrons are placed on the 
Cu-lattice sites |44j . 

Using the BCS gap equation for illustration, i.e. 



^ pp 2Ep' ■ 



(4) 



p' 



with Vppi , the effective interaction or Fourier-transformed vertex V pp 



= T(p' — p), we immediately see 
that the two pairing electrons avoid the repulsive parts of the Coulomb interaction by arranging their pair 



wave-function in a <L 



-orbit, with the simplest nearest-neighbor pairing given in momentum space by 



A p = A (cosp x — cospy). This is expressed by a kind of Pitaevskii-Brueckner relation, i.e. 



^ 2EL ' 



(5) 



This latter relation indeed seems to confirm PWA's conjecture in that the essence of the pairing mech- 
anism in low-T c SC is "dynamic screening", whereas in the high-T c cuprates it is "anisotropic momentum 
space pairing". In Fig. [2] we see - however - that the dynamics of the pairing interaction as reflected in 
the w-dependence of the gap function A(ui), is also of relevance for the high-T c cuprates. Here again, the 
real part of the gap function changes sign at a characteristic frequency <P related to spin fluctuations and 
at higher energies to p — d charge fluctuations. 



8 Will be inserted by the editor 

IV) Dynamics of the pairing interaction in the 1BH- and 3BH-models 



The question of whether one can speak of a "pairing glue" has recently been addressed by Maier, Poilblanc 
and Scalapino for the 1BH- and t-J models [30 . If the dynamics of the pairing interaction is due to "small- 
energy" (two-particle) excitations, such as the characteristic structures seen in the spin susceptibility, then 
one might speak of the interaction as being retarded on the relative time scales of interest and of a "spin- 
fluctuation" glue mediating d-wave pairing. This clearly is reminiscent of the usual phonon mediated 
pairing interaction in the low-T c SC [30J. 

It is well-known that for these latter systems the gap function <fi(k, u) or, more precisely, the anomalous 
part of the Nambu selfenergy 27(k;a;), is only very weakly depending on momentum k, corresponding to 
the local character of the pairing interaction. However, the dynamics of the gap function is important and 
it enters the dispersion (Cauchy) relations between the real and imaginary parts of </>(k, u>) — 0i(k, u>) + 
i0 2 (k,w), i.e. 

</>i(k,w) = - / — } dw (6) 

■K J U! — U! 

and, for w = 0, 



Mk,u, = 0)=2/M^0 da/ (7) 



A measure for the fractional contribution to the gap function </>(k, w = 0) and to pairing that comes 
from frequencies less than fi can then be defined [30 , i.e. 

It gives the relative contribution to the pairing of the retarded "glue" part and of the non-retarded, i.e. 
"instantaneous" part. In Ref [30_ it has been demonstrated that J(k, fi) is a useful diagnostic for the pairing 
glue (phonon contributions) in the case of Pb. /(k, fi) increases as fi passes through Wphonon + ^0; where 
Aq denotes the SC gap and u>o the characteristic Pb phonon frequencies. Its asymptotic value I(k, fi ^$> 
wphonon + ^o) exceeds unity reflecting the fact there exists an instantaneous Coulomb pseudo-potential. 
This leads to a negative, frequency-independent contribution (JJ'NR to 0(k, w) and at high frequencies 
J(k, fi) exceeds 1 by the "instantaneous" contribution — </>nr/</>(o)- Here, using numerical techniques, we 
examine the question of a "pairing glue" which offers a way of distinguishing different pairing mechanisms 
for the one-band and three-band Hubbard models in the relevant strong-correlation regime for the HTSC. 

The variational cluster approach (VCA) is particularly well suited for a study of the quantity I(k, fi). 
This is due to the fact that the VCA allows for accurately calculating the real and imaginary parts of the 
anomalous self-energy and, thus, of the gap function </>i(k, u>). This has already been shown in a recent 
letter, reproducing the experimentally found "gap dichotomy" of the nodal and anti-nodal gaps in HTSC 
as a function of doping [17] . 

In the VCA for the Hubbard model, the lattice is tiled up with (isolated) clusters of a given size. The 
corresponding Hamiltonian H' has the same on-site interaction U as the original Hubbard model, but it 
has modified single-particle hopping parameters t' (and chemical potential fi') with, in particular, if = 
between different clusters. The cluster provides a "reference system", which spans a space of trial self- 
energies E(H'). The self-energy that "best" describes the physics of the infinite lattice is then constructed 
via a variational principle, searching for the stationarity of the ground potential fi(S) in the subspace of 
the cluster self-energies S(H'). The trial self-energies S(H') are varied by varying the parameters t' , // 



Will be inserted by the editor 




Figure 4. (2x2) reference cluster used for 3BH-model. Each cell contains a CuC>2 unit. To reduce finite-size 
effects in 2-particle correlations, a more symmetric arrangement of O sites is implemented, compared to previous 
publications (see Ref. |25j). 



and the Weiss (i.e. pairing, AF, etc.) fields and a local chemical potential-shift term, which guarantees 
a consistent treatment of the particle density [23]. For the SC ground-state, a Nambu representation of 
the self-energy is used and 0(k, w) is obtained from the ri-componcnt. 

The results for J(k, fi) for the 3BH-model to be discussed below, are based on a (2x2) cluster with 4 
CuC>2 unit-cells and the self-energy (and Green's function) of this 12-site cluster was extracted from exact 
diagonalization. The geometry of our reference cluster is reproduced in Fig. |4j As we found 2-particle 
correlation functions on a small cluster to be more sensitive to symmetry than the 1-particle spectral 
functions and the phase diagram, we use a symmetric arrangement of the O sites in the reference cluster. 
Additionally, the Coulomb interaction U p d between holes on Cu and O sites is implemented with periodic 
boundary conditions. 

On the other hand, the results for the 3BH-model phase diagram in Figs. [5] and the 1-particle excita- 
tions in Figs.[6]a) & d) are in close agreement with our earlier results published in Ref. [25], where we used 
a (not-symmetrized) (2x2) cluster. The "upshot" of these results is, that the 3BH- and the 1BH- model 
phase diagram results with the latter presented in Figs, p^ a) & c) are very similar, reproducing in both 
cases the overall ground-state phase diagram of the high-T c superconductors (see also Fig. 3 of Ref. [25 ). 
In particular, they include salient features, such as the enhanced robustness of the AF state as a function 
of electron doping and the tendency towards phase separation ("PS" regime) into a mixed AF-SC phase 
at low doping and a pure SC phase at high (both hole and electron) doping. In the low-doping regimes, 
we find at a homogeneous symmetry-broken state in which both, the AF and the d-wave SC order pa- 
rameters m and A are nonzero. This corresponds to a phase "AF+SC", where AF and d-wave SC order 
microscopically coexist. A homogeneous state with pure d-wave SC (m = and A > 0) is obtained in the 
larger doping regimes. In the "in between doping" regions macroscopic phase separation occurs, where 
these two latter phases are thermodynamically unstable. In previous work, we have shown for the 1BH- 
model, that larger reference clusters still yield a qualitatively similar phase diagram, with one noteable 
exception: the larger clusters (up to 10 sites in the lBH-modcl) results suggest, that phase separation 
remains persistent in the h-doped case and disappears in the e-doped situation [28 . Based on the overall 
very similar phase diagrams for the (2x2) 3BH-model and for the (4x2) lBH-model (Figs.pjlb) & d)), 
we expect a similar disappearance of the "PS"-region for e-doping also in the 3BH-model. 



Why are the phase diagrams in Fig. [5] so similar despite the fact that, for example, the p-d charge 
fluctuations give rise to an at higher energies (w ~ Z\ pd ) even qualitatively different behavior of the 
dynamic pairing interaction for the 3BH results when compared to the lBH-data as shown in Figs. [6] & [7] 
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Figure 5. Ground-state phase diagram for the 3BH- model ( a) & b) ), based on a (2x2) reference cluster with the 
parameters specified in the text, in comparison with the lBH-model ( c) - f ) ), based on a (4x2) reference cluster. 
In d), e) & f), the chemical potential \i is plotted as a function of doping. The corresponding order parameters for 
the models are given in ( a), c) & e) ). Between the mixed phase with both antiferromagnetic and superconducting 
order (AF+SC) and the pure superconducting phase (SC) at higher-dopings, is a region of phase separation (PS), 
described in text. For lBH-model, both hole- and electron-doping is presented. There is no qualitative difference 
between 1BH and 3BH phase diagrams. 



for /(k, J?)? The reason is that on the scale of the (maximal) d-wave gap energies, A p< i is a doping- 
independent "high-energy" scale (~ ten times larger than the SC gap) the role of which is taken over 
in the 1BH- model by another doping independent "high-energy" scale, namely an effective Hubbard U 
(see discussions, below). As just stressed, however, the dynamics of the pairing mechanism, may be 
quite different, which is what we find indeed below. We expect this to be of importance for the "non- 
universality", i.e. material-dependence of the cuprate SC 



Also the single-particle excitations (as displayed in Figs. 6 up to 7 as well as in our earlier Ref. 
[25J) are found qualitatively to be similar, concerning the "low-energy" physics and, in particular, the 
e-h asymmetry. This asymmetry explains the corresponding e-h asymmetry found in the ground-state 
phase diagram for the robustness of the AF phase (Figs. [5] c) & e): doped holes first enter around the 
nodal point (?, %), where the SC gap vanishes, introducing a "gap-less" screening which very effectively 
destroys long-range AF order [TB]. In contrast, introducing electrons around the anti-nodal point (tt,0) 
fixes jx (due to the large density of states) within the SC gap, in a regime where this gap is maximal. 
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Figure 6. 3BH-model: single-particle spectral function A(k,cu) ( a) & d) ), d-wave projection Kj(w) of the spin 
susceptibility ( b) & e), see eq. (10)) and integrated contribution of gap function J(k = (-7T, O.I71"), Q) ( c) & f) ) 
for a (2x2) reference cluster (as in Fig. 4). For a) - c), the hole-doping is 4% and for d) - f), the hole-doping is 
15%. For both dopings, the system is in a pure superconducting state. 
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Figure 7. lBH-model: single-particle spectral function A(k, w) ( a) ), d-wave projection Vd(u>) of the spin 
susceptibility ( b), see eq. (10) ) and integrated contribution of gap function J(k = (7r,0.27r), J?) ( c) ) with a 
(3x3) reference cluster. Hole-doping is 15% and the system is in a pure superconducting state (see Fig. pi. 



Here, an incomplete screening cannot disrupt AF order up to significantly large e-dopings where, finally, /i 
enters also here the "gap-less" (f , f ) region. So again, it is the corresponding low-energy physics embedded 
in the qualitatively rather similar single-particle spectra plotted in Figs. 6-7, which determines salient 
(here AF) features of the phase diagram. 

So, let us finally come back to the question of the dynamics in the pairing interaction and the fractional 
contribution /(k, I?) to the gap function in the 3BH- as well as 1BH- models. 

Our results for /(k, ft) for the 3BH-model, i.e. for the fraction of the zero-frequency gap function 
which arises from frequencies below i?, are plotted in Figs. [6]c) & f). Here two doping cases, 4% hole 
doping and 15% hole doping are shown, together with the corresponding single-particle spectral function 
A(k,uj) and the data for the d-wave projection of the cluster spin susceptibility (eq. (10)). 



J(k, Q) is plotted for the 3BH-model for a doping x = 4% in Fig. [61c) and for a higher doping x — 15% 
in Fig. [6]f). k = 7r(0.1,0), i. e. a k-point close to the anti-nodal point (n,0) where the d-wave SC gap is 
maximal. For both dopings, we observe a first steep rise at a typical energy <P of about uj + Z\ , where 
Aq denotes the quasi-particle (SC) gap and ujq a characteristic spin-fluctuation frequency. Close to half- 
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Figure 8. Integrated contribution of gap function J(k = (n, 0.l7r), fi) for the 3BH model. In a), the first dominant 
increase in I(k, fi) is plotted for different values of Udd/t. In b), the doping dependence of 7(k, fi) is shown, with 
the modest increase of p-d charge fluctuations as function of doping. 



filling, ujq is roughly estimated by the strong-coupling result for the Cu-Cu exchange energy luq = 2Jc Ul 
with 



Jcu — 



Kd 



(A pd + U pd ) 2 A pd + Upp/2 U dd 



= 0.lt 



pd • 



(9) 



As shown by the d-wave projected result Vd,d(w) in Fig. m\ b) & e), for the imaginary part of the 
cluster spin susceptibility, with the definition 



1 






V d ,ci(u) = -^2 ^_j (cos(fc ;c ) - cos(fe v )) (cos(%) - cos(fc^)) Imxd(k- k',w), 

k,k' 



(10) 



the first dominant peak is found at low doping, i.e. 4% at about this energy (here ujq 
in units of the hopping t pd ) and, at 15% doping at a higher energy ujq = 0.2. 



0.13 as always 



7(k, J?) displays an amazingly prominent first increase at an energy <S> of about Aq + ojq, where Aq 
is the quasiparticle gap and loq the first dominant peak in the d-wave projection of the dynamic spin 
susceptibility Im/j^k, w). The quasiparticle gap is found to be at A = 0.52 (obtained from the spectral 
function) for 4% doping and at a similar value (Z\o = 0.41) for 15% doping, thus Aq + ljq = 0.65. The 
first sharp rise in 7(k, Q) at the low-doping situation (~ 4%) is at about <P = 0.56. We observe that in 
this low-doping case for the 3BH-model more than 90% of the saturation value of I(k, fi) = 1 is due to 
the dynamic contribution of the spin fluctuations. 

For 15% doping ojq is increasing [13] , and so is wo + ^io- Fig. ^a) further confirms the leading role of 
the spin-fluctuations. Here the first dominant increase in 7(k, fi)is plotted as a function of U/t at 15% 
h-doping. One sees that as the value of U/t increases, this dominant increase in J(k, fi) is shifted to lower 
energies, scaling essentially with the exchange energy J(Aq is found to be essentially constant as function 
of U/t) Here, at 15% doping, the spin-fluctuation contribution to the saturation value of 7(k, Q) = 1 
slightly decreases compared to the 4%-doping situation. From there on a mild rise takes place, starting at 
about u = 1 for both dopings. In both doping cases this corresponds to virtual e-h transitions involving 
the charge transfer gap ~ A pd /2 (see Fig. [7|i 



Fig. [8] b) illustrates the influence of the doping dependence of p-d charge fluctuations on 7(k, ]?). 
Here, at first glance, the influence of these fluctuations on the integrated gap function 7(k, fi) appears 
to be very mild. 
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However, this picture changes significantly, when we confront the 3BH-model results for I(k, fi) in 
Figs. [6]d) - f) with the corresponding ones for the lBH-model and, thus, the doped Mott-insulator case 
(Figs!"[7m - c): again, we find for the latter model a first dominant rise in the integrated gap function. (We 
plot here only the results for the 15%-doping case with the low-doping (4%) results being qualitatively 
similar.) This first rise - when analyzing its ([//independence - can again be attributed to dynamic 
spin excitations. However, a second prominent increase in this lBH-model case occurs after a rather 
pronounced drop at energies of order (U/t)/2. This has already been observed in this lBH-model case in 
Ref. [30J and was termed a "non-retarded" contribution occuring at an energy scale set by the Mott gap 
and being related to excited states (see the red arrow in the corresponding A(k,w) spectrum in Fig. [Ha) 
involving the upper Hubbard band. For the t-J model, Ref [3U] finds this energy scale being pushed 
to "infinity", with the corresponding exchange interaction being instantaneous. In contrast, the retarded 
(spin-fluctuation) interaction , i.e. the first steep rise in /(k, fi), occurs at an energy scale which is small 
(cjo — O.li) compared to the bare bandwidth 8i (and the bare U). Thus, the lBH-model as well as the 
t-J model may be interpreted to contain both retarded and non-retarded contributions. 

This is qualitatively different in the 3BH-model: Here, in Figs.[6]c) and f) and Fig.[8]b) we see a rather 
continuous "filling in" of integrated weight in the Q regions in I(k, J?), where in the lBH-model (Fig.[7]c) 
) this weight is clearly missing. This "filling in" is due to electron-hole excitations of 0(2t p d) and, thus, 
due to p-d charge fluctuations (see the corresponding "red arrow" for the e-h transitions in A(k,uj)). 

Thus, in the sense of the interpretations used for the lBH-model above (and in Ref J30J) we may 
conclude that the 3BH-model is indeed different and contains only retarded contributions in its d-wave 
pairing interaction. 

In summary, then, the question posed in sect. Ill, i.e. whether there is a "pairing glue" in the 3BH- 
model is a question of whether the dominant contribution to the pairing function </>i(k, w = 0) comes 
from the integral of <p2(k,u>)/uj and, more specifically, a "low-energy" (compared to 8t) region in this 
integral. Earlier results have shown that both the 1BH- and t-J models exhibit spin-fluctuation "pairing 
glue" [30 31j. For both the 3BH-model as well as for the lBH-model results, with the latter displayed for a 
similar higher-doping case (15%) in Fig. [7] one can include the full range of virtual e-h transitions between 
the lower and upper Hubbard bands. As a consequence the dispersion (Cauchy) relation for 0(k, uj = 0) 
does, strictly speaking, not contain a non-retarded contribution </>nr and the asymptotic value of I(f2) 
approaches unity. In the lBH-model, we observe a first steep rise, which can again be identified with a 
typical spin-fluctuation energy. Then, at higher energies after a drop a second steep rise occurs at energies 
corresponding to (virtual)e-h interband transition of order ~ ^. As argued in Ref. [3D], it is this part 
which in the t-J model for U — > oo (i. e. when the upper Hubbard band is projected out) gives rise to 
an instantaneous contribution </>nr- Thus, in this sense, both 1-band models also display a non-retarded 
interaction. 

Our new result here is that in the SBH-model only low-energy (spin-fluctuation) retarded contributions 
to the pairing interaction dominate, whereas in the lBH-model also a significant high-energy contribution 
occurs. While these differences between a Mott and a charge-transfer insulator appear to be renormalized 
away in the very similar, i.e. "universal" ground-state phase diagram of these two models, one may use 
Eliashberg-type of arguments to form an expectation that the finite-T phase is rather different in the two 
models, i.e. strongly material-dependent (see also Ref. [27 J. 



V) Spin susceptibility 



As is clear from the previous sections, two-particle (2-p) excitations, especially magnetic ones, are fun- 
damental for understanding the issue of the "pairing glue" in high-T c cuprate superconductors. In addi- 
tion, for hole-doped cuprates in the superconducting (SC) state, the celebrated resonant magnetic mode 
emerges with its peak intensity being highest around the wave vector Qaf — (tt, tt) |45j . Away from Qaf, 
the mode has both a downward and upward "hour-glass"-like dispersion. For electron-doped materials |4l)] 
the magnetic excitation spectrum has a different structure: it is confined to a small momentum region 
around (jr, 7r) and it is essentially dispersionless |35| . 
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In order to provide an appropriate description of the magnetic susceptibility and, in particular, of the 
magnetic resonance mode it is important to adopt a theory which is working in the experimentally relevant 
strong-correlation regime. As discussed above, VCA provides a suitable tool to study strongly-correlated 
systems for large lattices. In Ref. (21] we have developed a novel theory for two-particle excitations based 
on the VCA. This approach is, in contrast to previous weak-coupling and/or semiphenomenological 
treatments, parameter free and applies to the relevant strong-coupling regime of the one- and three-band 
Hubbard models. Within this method, two-particle susceptibilities are evaluated by first computing the 
exact two-particle irreducible vertex r within the reference system (cluster), and then inserting T into 
the Bethe-Salpeter (BS) equation for the exact (lattice) two-particle susceptibility [M]. Notice that the 
BS equation is in principle exact and not restricted to weak-coupling (in contrast to RPA). In principle, 
r depends on three frequencies and four cluster indices. Therefore, solving the BS equation with the 
corresponding full frequency and momentum dependence of i" 1 is a numerically quite expensive task. For 
this reason, we replace the full _T with its average over external frequencies, and set the site indices pairwise 
equal, which corresponds to averaging over the external cluster momenta. After this simplification, r(iuj) 
can be computed as [34147] 

r(i w )=[(x°'M)- 1 -(x , M)- 1 ] , (ii) 

where x' is the exact and x° is the "bubble" susceptibility of the reference system, and w is a Matsubara 
frequency. The term x' in eq. (11) is evaluated directly within a Lanczos procedure and the first one 
is computed using the exact cluster Green's functions. The effect of the above approximations is then 
partially compensated by multiplying r(Q 1 ico) by a constant term a which is determined so as to fulfill 
the sum rule for the transverse spin susceptibility. Details can be found in Ref. [34 . Finally, the lattice 
spin (or, corresponding charge) susceptibility, obtained by solving the BS equation, is given by 

X (Q, ioj))- 1 = X °(Q, iw)) _1 - aT(iu) , (12) 

where X°(Q:* W )) _1 is the lattice "bubble" susceptibility evaluated with the fully dressed lattice Green's 
functions. Since we are studying the superconducting phase, the bubble susceptibility has a contribution 
from both normal and anomalous Green's functions. |34| Due to the partial breaking of translation 
symmetry introduced by the cluster tiling, both the xo an d X depend on the wave vector Q of the 
correspondingly reduced Brillouin zone, and are matrices in Nambu space as well as in the cluster sites. 

In Ref. [25 , as well as in the present paper, the qualitative similarity of the phase diagram and single- 
particle excitations between the three-band and the single-band Hubbard models was demonstrated. In 
particular, it was explicitly confirmed that the asymmetry between electron- and hole-doped cuprates, 
despite being of fundamentally different nature, shows very similar signatures in the single-particle spec- 
trum of both models, provided a next-nearest-neighbor hopping t' , is included in the latter. The question 
remains whether this similarity can be extended to two-particle excitations. To address this question, 
we compare in Fig. [9] the imaginary parts of the spin susceptibility evaluated in the deeply underdoped 
regime [x = 4%), for both the single- and the three-band models. 

While there are some minor differences, both spectra display a dispersion which is a remnant of 
the spin- wave band in the antiferromagnetic half-filled phase. This fact signals the presence of strong 
antiferromagnetic fluctuations at this doping concentration also in a doped charge-transfer insulator. 

Last not least, the VCA- extension to 2-particle susceptibilities can account for - without any adjustable 
parameters - the magnetic resonance in both hole- and electron-doping regimes. Close to optimal doping, 



the magnetic spectrum for the h-doped single-band Hubbard model (Fig. 10 a) ) shows the famous 
magnetic resonance peak at energies of about O.li — 0.2i and around wave-vector (jt, it). From this 
resonance, a downward dispersion extends down to the onset of the particle-hole continuum. The striking 
feature leading to the celebrated "hourglass" structure is, however, the additional "upward" dispersion, 



which can be seen in the spectrum of the single-band Hubbard model displayed in Fig. 10 a). Because of 
computational difficulties, we were not able to obtain this feature in the three-band model calculation. 
This is due to the fact that the smaller (2x2) CuC^-units cluster adopted for the three-band model still 
has significant antiferromagnetic correlations which prevent the formation of the hourglass structure (a 



similar observation holds also for the lBH-model, where we had to go to (3x3) clusters (Fig. 10) to be 
able to observe the resonance). 



In the electron-doped case, the magnetic excitation spectrum displays a different structure. In Fig. 10 b), 
we present the first strong-coupling and parameter-free evaluation of the magnetic spectrum in the 
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Figure 9. Imaginary part of the spin susceptibility \ s , w) f° r the single-band ( a) ) and for the three-band ( 
b) ) Hubbard model at doping x = 4%. Results are obtained with a reference system of (2x2) unit cells. 
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Figure 10. Imaginary part of the spin susceptibility for the single-band Hubbard model ((3x3) reference clusters, 
U — 8t) in the hole-doped (x — 0.15, a) , taken from Ref. 34_ ) and in the electron-doped (x — 0.14, b) ) cases. 



electron-doped case, within our VCA approach for two-particle excitations. As one can see, in this case 
the weight of the magnetic spectrum is concentrated essentially in the vicinity of (7r,7r), in accordance 
with the experimental situation. The downward dispersion observed in the calculation only carries very 
little weight, when compared with the strong peak around (tt,tt). Therefore, this might be the reason 
why this dispersion is not seen in inelastic neutron scattering experiments [46 . Finally, in contrast to the 
hole-doped case, the upward dispersion is not observed for electron doping. This is also in accordance 
with the experimental observations |35| . 

We have also extended, in a similar spirit, the dynamical cluster approximation (DCA) to 2-p suscep- 
tibility calculations at finite-temperatures [17] for the infinite-lattice case. In Ref. [UJ, this scheme has 
been used to obtain - in the case of the lBH-model - information on the finite-T phase diagram (via the 
instability towards a symmetry-broken phase of the corresponding susceptibility). 



VI) Conclusion 



One key issue in the research on high-T c cuprate superconductivity is to derive a detailed understanding 
for the specific role of the Cu-d x i_ y i and 0-p x , y orbital degrees of freedom. Very early on a corresponding 
three-band Hubbard (3BH) model has been suggested but in most studies an effective one-band Hubbard 
(1BH) or t-J model has been used. The one-band models have been shown, using a variety of techniques, 
to reproduce salient features of the cuprates such as the competing antiferromagnetic (AF) and super- 
conducting (SC) phases including the electron-versus-hole doping asymmetries. Here, the physics is that 
of doping into a Mott-Hubbard insulator, whereas the actual high-T c cuprates are doped charge-transfer 
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insulators. Using a cluster embedding scheme, i.e. the variational cluster approach (VCA) which allows 
for a controlled route to go to the large system size limit (and, thereby, low-energy limit), we compare in 
detail the competing phases in the ground-state, the single-particle and two-particle excitations and, in 
particular, the dynamic pairing interaction and the issue of a "pairing glue" of the 3BH and 1BH model. 
This study demonstrates that there are pronounced differences on a "higher-energy" scale (compared 
to the SC gap energy) for the two models, such as the charge-transfer energy A p d in the 3BH-model. 
While these "high-energy" features, where the material-dependent physics enters, play no decisive role for 
the overall structure of the phase diagram, which shows excellent agreement between 3BH- and 1BH- 
models, there appear significant differences in the dynamics of the pairing interaction: in the 3BH-modcl 
the interaction is dominated by a retarded pairing due to low-energy spin fluctuations, whereas in the 
1BH- model in addition a part comes from excited states involving the upper Hubbard band, which may 
be associated with a "non-retarded" contribution. Thus, one may term the ground-state phase diagram 
"universal", whereas the different dynamics of the pairing interaction of the two models should be reflected 
in a pronounced material dependence of the finite-T phase diagram. 

First steps in this direction have already been undertaken in Ref. |27] and are also presently followed 
up in our group |48| . 

Finally a comment about the RVB-type of physics in the two models: A v ^ also enters in the exchange 
coupling J between spins on the Cu sites. Thus, the role of U in the 1-band model, setting the char- 
acteristic energy of the RVB-coupling J, is taken over (see Eq. (9)) in the 3-band model by the energy 
of the p — d charge fluctuations. In that sense, one may argue, that the final increase in the relative 
pairing strength /(k, w) to the value of I(k,cu) = 1 again reflects RVB physics. However, our results 
demonstrate that there is quite a difference with respect to the relative weights of the spin-fluctuation 
and higher-energy contributions between the two models. 
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